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Using a circuit QED device, we demonstrate a simple qubit measurement pulse shape that yields 
fast ring-up and ring-down of the readout resonator regardless of the qubit state. The pulse differs 
from a square pulse only by the inclusion of additional constant-amplitude segments designed to 
effect a rapid transition from one steady-state population to another. Using a Ramsey experiment 
performed shortly after the measurement pulse to quantify the residual population, we find that 
compared to a square pulse followed by a delay, this pulse shape reduces the timescale for cavity 
ring-down by more than twice the cavity time constant. At low drive powers, this performance is 
achieved using pulse parameters calculated from a linear cavity model; at higher powers, empirical 
optimization of the pulse parameters leads to similar performance. 


Over the last decade, circuit quantum electrodynam¬ 
ics [I] (cQED) has become a leading architecture for con¬ 
structing scalable networks of solid-state qubits, finding 
application in the context of not only superconducting 
qubits 0J but also spin qubits d] and potentially other 
systems [5|. In this paradigm, each qubit is coupled to a 
resonator in which it induces a state-dependent frequency 
shift, allowing the qubit state to be interrogated using a 
pulsed tone near the resonator frequency. A great deal of 
research has focused on optimizing the speed and fidelity 
of such measurements. Most significantly, the ongoing 
development of quantum-limited amplifiers has im¬ 
proved achievable signal-to-noise ratios enormously. The 
introduction of Purcell filters [3 El lias enabled the 
use of resonators with fast time constants, whose high 
bandwidth would otherwise provide a pathway for qubit 
relaxation vi a sp ontaneous emission (Purcell effect Il2|). 
Some work 11 M has also explored the use of pulse 
shapes with an initial overshoot in order to populate the 
readout resonator more quickly than the standard square 
pulse. 

Although these improvements have made fast, high- 
fidelity qubit readout in cQED systems routine, rela¬ 
tively little attention has been devoted to the problem 
of returning the readout resonator to its ground state 
immediately after the measurement pulse. If the pulse 
is simply turned off, residual photons gradually exiting 
the resonator will continue to measure the qubit 14], 
preventing high-fidelity operations for a period of several 
time constants. Even for a resonator with a fast time con¬ 
stant, this delay is typically longer than the time needed 
for qubit control and measurement. A technique for re¬ 
ducing the residual population on a timescale faster than 
the resonator’s free decay is therefore desirable in any al¬ 
gorithm in which qubits need to be re-used shortly after 
measurement, e.g. error correction with surface |15lil6|. 
C4 171 or Bacon-Shor El codes. A major impediment 
has been the fact any such scheme needs to work in the 
absence of prior knowledge as to which of the two possible 
state-dependent resonator frequencies will be realized. 

Here we present the first experimental demonstration 


of driven state-independent reset of a readout resonator, 
using a specially designed yet simple pulse shape that we 
term the Cavity Level Excitation and Reset (CLEAR) 
pulse [Fig. Oja)]. The pulse uses short segments to “kick” 
the resonator rapidly from one steady-state population to 
another: at the beginning of the pulse, two such segments 
drive the population from zero to the desired steady-state 
value, and at the end, two more drive it back to zero. In 
this work, we focus on quantifying the effectiveness of 
the depopulating segments. Using a Ramsey experiment 
to extract the number of residual photons in the cavity 
after the pulse, we compare the performance to that of 
a standard square pulse. We find that for pulse powers 
where the cavity response remains linear, the theoreti¬ 
cally derived CLEAR pulse shape depopulates the cavity 
to a negligible level in a time more than two cavity time 
constants faster than that needed after the square pulse. 
At higher powers, optimizing the pulse shape empirically 
using an iterative algorithm leads to equally good perfor¬ 
mance. 

The experimental device is a fixed-frequency trans- 
mon qubit mounted in a 3D aluminum cavity ITbi at¬ 
tached to the mixing chamber of a dilution refrigera¬ 
tor at an indicated base temperature of 10 mK. Qubit 
and measurement drive tones are generated using Ag¬ 
ilent E8267D function generators and modulated using 
Tektronix AWG7000 series arbitrary waveform genera¬ 
tors at a 2 GS/s sample rate. Qubit pulses are 4 a Gaus- 
sians with DRAG [20j] correction. The cavity is measured 
in transmission, and the transmitted signal is fed to a 
HEMT amplifier (Low Noise Factory LNF-LNC6_20A) 
at 4 K using a superconducting NbTi/NbTi semi-rigid 
coaxial cable. After additional amplification at room 
temperature, the signal is mixed down to 16 MHz and 
digitally demodulated. The cavity is measured to have 
bare frequency /bare = 10.7457 GHz, dressed frequency 
/dressed = 10.7594 GHz, and linewidth n/2n = 1.1 MHz 
(corresponding to a time constant T cav = 1/k = 0.14 ^s). 
The qubit has frequency /oi = 4.83315 GHz, anhar- 
monicity 5/2n = —155 MHz, average Xj « 50 \x s, and 
average Tf cho ss 60 /rs. Preparing the qubit in the excited 
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state shifts /dressed by the cavity pull 2x/2n = —2.6 MHz; 
the measurement tone is applied at the midpoint of the 
two frequencies, /dressed - x/ 2 ^- 

The residual population after a measurement pulse can 
be quantified in terms of the mean cavity photon number 
n at some time after the end of the pulse. We use the 
sequence illustrated in Fig. QIb) to extract n following 
an initial measurement pulse denoted Ml. The measure¬ 
ment pulse is followed by a quick Ramsey experiment 
(/gate = 8 ns, tR = 0 to 600 ns) to probe the ac Stark 
shift and dephasing from any residual photons 14|, 21]. 
The time R e i ax between the end of Ml and the start of 
the Ramsey experiment can be varied to measure n as 
a function of time after the end of Ml. The time fbuffer 
between the Ramsey experiment and the second mea¬ 
surement pulse (M2) is set to 400 ns to ensure that even 
when tji and R e i ax are both short, M2 is not corrupted 
by lingering photons from Ml. 


A typical Ramsey trace is shown in Fig.[2]]a). The non¬ 
monotonic modulation in both amplitude and frequency 
arises from the fact that the cavity population evolves 
during the Ramsey delay, leading to recurrences. We 
derive the expected form of the Ramsey trace during this 
transient response usin g th e positive-R-function method 
as in Gambetta et al. [14], where it was applied to the 
steady-state problem. For a Ramsey detuning A (here 
10 MHz), decoherence rate T 2 , and initial phase </>o, the 



FIG. 1. (color) (a) Schematic shape of the piecewise-constant 
CLEAR pulse (solid red line) with a square pulse (dashed 
black line) for reference, (b) Pulse sequence used to extract 
the residual cavity population after a measurement pulse. The 
qubit is prepared in either the ground or excited state and 
then a measurement pulse (either a square pulse or a CLEAR 
pulse) of length tMi is applied. After the measurement pulse, 
an adjustable delay Relax precedes a pair of Xqq pulses sepa¬ 
rated by tn comprising a Ramsey experiment. The Ramsey 
experiment is followed by another brief delay tbuffer and a 
square measurement pulse of length tM 2 = 10 ps. 


resulting functional form is 


S{t R ) = ^[1 - Im(exp(—(r 2 + A i)t R + (</> 0 - 2n 0 xr)i))], 

(!) 

where t = (1 — e ( K+2x ^ tR ) / (n + 2\i) and no is the value 
of n at the beginning of the Ramsey experiment. Using 
k and x obtained from frequency-domain measurements 
and taking T 2 = l/T| cho , the only free parameters are no 
and <j> o- As illustrated in the figure, this function yields 
a good fit to the data, allowing reliable determination of 
no. In the rest of this work, we use the extracted no to 
quantify the residual population as a function of pulse 
shape, drive power, and wait time. We note that no 
does not include the background thermal population of 
the cavity, which is accounted for by the T 2 term and is 
calculated from Tf cho to be ~ 0.02 on average (assuming 
thermal photons are the dominant source of steady-state 
dephasing [HSl). 

The Ramsey fit method of obtaining no was validated 
by using it to measure no as a function of both wait 
time Relax and pulse power. For simplicity, these tests 
were performed using square pulses for both Ml and M2, 
varying only the power of Ml, denoted P as illustrated 



FIG. 2. (color) (a) Sample Ramsey experiment and fit. The 
black curve is a fit of Eq. [l]to the data (red circles), yielding 
initial cavity population no ss 0.9. (b) Markers indicate no 
versus wait time Relax after a square measurement pulse with 
drive power P norm = 2. Here and throughout, blue and red 
denote experiments in which the qubit was prepared in the 
ground and excited states, respectively. Solid blue and red 
curves are exponential fits to the respective data sets, (c) no 
versus drive power measured at Relax = 40 ns after a square 
measurement pulse. The dashed line indicates the prediction 
for a linear cavity, accounting for Relax- The solid curves 
account for the cavity’s self-Kerr nonlinearity calculated from 
the measured qubit and cavity parameters. 
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in Fig. [Ha). For convenience we define the normalized 
drive power P norm = P/Piph, where Pi p h is the steady- 
state drive power that yields n = 1 , as inferred from a 
standard Ramsey experiment measuring the Stark shift 
A u) = 2 xn 14] induced by a CW tone. Figure [H[b) shows 
no extracted from Ramsey fits as a function of t re iax- Re¬ 
gardless of the prepared qubit state, the decay is expo¬ 
nential, as expected for free decay, and the time constant 
T cav extracted from the best-fit curve is consistent with 
the value of n obtained from frequency-domain measure¬ 
ments. Figure [H]c) shows no as a function of P norm at 
treiax = 40 ns; similar behavior was observed at all values 
of Peiax for which non-negligble no were measured. The 
dashed line indicates the expected behavior assuming a 
linear cavity: no = e _Ktrelax P norm . The data exhibit a 
transition from a linear response at low powers to a super- 
(sub-) linear response at high powers when the qubit is 
prepared in the ground (excited) state. This behavior is 
consistent with the cavity’s expected self-Kerr nonlinear¬ 
ity 0, which shifts the cavity frequency by a negative 
amount K per cavity photon, pushing it closer to (far¬ 
ther from) the measurement frequency when the qubit is 
in the ground (excited) state. Approximating K in the 
small-<5 limit as I\ = 2g 4 d(3u 4 + 2uj q io^ + 3Lo 4 )/(uj q —ij^) 4 , 
where uj q = 27t/oi, uj r = 27r/d reS sed, and g is calculated 
as in Ref. [Hj], we obtain K ss — 14 kHz. The solid blue 
(red) curve indicates the expected cavity response with 
the qubit in the ground (excited) state, accounting for 
the calculated non-linearity. 

Having thus validated our method of quantifying resid¬ 
ual photons, we then switched to a CLEAR pulse shape 
for All; for consistency, we continued to use a square 
pulse for M2. For an ideal single qubit-cavity system, 
the optimal CLEAR pulse envelope [Fig. Oja)] consists of 
five piecewise-constant segments: two ring-up segments, 
one steady-state segment, and two ring-down segments. 
For the pulse to behave as intended, its bandwidth needs 
only to be much greater than that of the cavity, a con¬ 
dition readily achieved in our setup. Setting the carrier 
frequency to /dressed ~ x/27t, as done here, is not required 
but maximizes both SNR and simplicity: in this case, a 
given measurement pulse yields the same steady-state n 
regardless of qubit state as long as the cavity remains in 
the linear-response regime. The lengths of the ring-up 
and ring-down segments were all initially fixed at 150 ns 
(approximately T cav ), and their amplitudes relative to 
that of the steady-state segment were calculated by solv¬ 
ing a driven damped harmonic oscillator model to find 
the pulse shape that populates and depopulates the cav¬ 
ity in the shortest amount of time regardless of the qubit 
state. 

The cavity IQ plane trajectories produced by square 
and CLEAR pulses with the same steady-state ampli¬ 
tude (Pnorm = 3.6) are shown by markers in Figs. [3](a,b) , 
respectively. For all trajectories, the time step between 
markers is 24 ns. Solid lines in these plots indicate the 


theoretically calculated response of the cavity to each 
pulse, multiplied by an overall amplitude factor to match 
the data and adjusted to reflect the independently ob¬ 
served 20% thermal population of the qubit excited state 
(which was reduced by an order of magnitude on later 
cooldowns of this device with additional input line atten¬ 
uation). We see that the experimental cavity responses 
track the theoretically calculated ones very well, and that 
compared to the square pulse, the CLEAR pulse yields 
more compact trajectories that reach near-steady-state 
populations (at both n ~ 3.6 and n ~ 0) in less time. 

We quantitatively compare the performance of the 
CLEAR pulse to that of a square pulse using the Ram¬ 
sey fit method. To provide a fair comparison, a zero- 
amplitude segment is appended to the square pulse to 
allow undriven decay during a time equivalent to the to¬ 
tal length of the CLEAR pulse’s two ring-down segments. 
The results are shown in Fig. He). At all measurement 



FIG. 3. (color) (a,b) IQ-plane cavity trajectories in response 
to a square pulse and CLEAR pulse, respectively. In each 
plot, experimental results (markers) are superimposed on the¬ 
oretical calculations (solid curves), the dashed circle indicates 
the target population n = 3.6, the black cross indicates the 
origin, and the black arrows indicate the directions of the tra¬ 
jectories. (c) Residual cavity population versus pulse power 
for both square and CLEAR pulse shapes. For the CLEAR 
pulse, the Ramsey experiment begins immediately at the end 
of the pulse, while for the square pulse, a delay of approx¬ 
imately 300 ns is inserted to match the total length of the 
CLEAR pulse’s two ring-down segments. 
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powers, the CLEAR pulse significantly outperforms the 
square pulse; moreover, for drive powers that keep the 
cavity in the linear regime (evidenced by no being inde¬ 
pendent of the prepared qubit state), the residual popula¬ 
tion immediately after the CLEAR pulse is negligible. As 
seen in Fig. 0(b), for drive powers in this range, negligible 
no is not obtained until ~ 600 ns after a square pulse; al¬ 
lowing for the 300 ns taken by the ring-down segments of 
the CLEAR pulse, we find a net speedup of ~ 300 ns, or 
approximately 2T cav . At higher powers, it appears that 
the cavity nonlinearity, not taken into account in calcu¬ 
lating the optimal CLEAR pulse parameters, prevents 
perfect ring-down. We also find non-idealities when the 
lengths of the ring-down segments are reduced in an effort 
to shorten the ring-down time: for 120 ns ring-down seg¬ 
ments (a 20% reduction), we find measurable no even in 
the linear regime, as illustrated in Fig. |U(a). Further re¬ 
ductions in the segment lengths increase the performance 
degradation. 
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FIG. 4. (color) (a) no versus drive power for the shortened 
CLEAR pulse (120 ns ring-down segments), for t re i ax = 0. (b) 
Evolution of no at each step of an empirical optimization algo¬ 
rithm for the shortened CLEAR pulse with P n orm = 10. Inset: 
shortened CLEAR pulse shape before optimization (magenta) 
and after (green), (c) Ramsey traces obtained using initial 
shortened CLEAR pulse, yielding no ~ 2.2 for ground and 
no ~ 0.91 for excited, (d) Ramsey traces obtained using fi¬ 
nal shortened CLEAR pulse, yielding no < 0.1 for both qubit 
states. 

To improve performance of the CLEAR pulse both at 
high powers and with shortened ringdown segments, we 
use an empirical technique to optimize the pulse parame¬ 
ters. As an example, we take as a starting point the pulse 
with 120 ns ring-down segments and a steady-state drive 
power of Pnorm = 10 (which yielded no ~ 2.16 for the 
ground state and no ~ 0.91 for the excited state), and 
run an iterative optimization algorithm that attempts 


to minimize no by adjusting the amplitudes of the ring- 
down segments. We keep f re iax = 0 throughout. The 
evolution of no with each iteration is shown in Fig. SKb): 
in fewer than 300 iterations, the pulse is optimized to 
yield no < 0.1 regardless of initial qubit state. Ramsey 
experiments before and after running this optimization 
are shown in Fig. [2c, d), revealing that coherence is pre¬ 
served with the optimized pulse shape regardless of initial 
qubit state. As seen in the inset of Fig. |U(b), the opti¬ 
mization process significantly increases the amplitudes of 
both ring-down segments. Extending our theoretical cal¬ 
culations to the non-linear regime may shed light on this 
result and potentially eliminate the need for empirical 
tune-up in this regime. 


In summary, using a single-qubit cQED system, a 
qubit-state-independent reduction in the time needed 
to reach a steady-state resonator population both dur¬ 
ing and after a qubit measurement pulse was achieved 
by including extra constant-amplitude segments in the 
pulse. For low-power drives, near-perfect ring-down 
(quantified using Ramsey experiments) is achieved us¬ 
ing segment amplitudes calculated from system param¬ 
eters. At higher drive amplitudes, similar performance 
is obtained following empirical optimization of the pulse 
shape. Though this demonstration used a 3D trans- 
mon, the same technique should be applicable to any 
cQED system; it may also be combined with machine¬ 
learning based analysis [25] and Purcell filters to fur¬ 
ther reduce the measurement cycle time. Future areas 
of interest may include numerical calculation of the opti¬ 
mal CLEAR parameters in the non-linear regime, exten¬ 
sion of this technique to resonators coupled to multiple 
qubits, and investigation as a possible method for im¬ 
plementing the resonator-induced phase (RIP) gate j26j 
non-adiabatically. 
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